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Laser driven plasma accelerators promise much shorter particle accelerators but their development 
requires detailed simulations that challenge or exceed current capabilities. We report the first direct 
simulations of stages up to 1 TeV from simulations using a Lorentz boosted calculation frame 
resulting in a million times speedup, thanks to a frame boost as high as 7 = 1300. Effects of the 
hyperbolic rotation in Minkowski space resulting from the frame boost on the laser propagation in 
the plasma is shown to be key in the mitigation of a numerical instability that was limiting previous 
attempts. 
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Laser driven plasma waves produce accelerating gra- 
dients orders of magnitude greater than standard accel- 
erating structures (which are limited by electrical break- 
down) [TJ [2] . High quality electron beams of energy up-to 
1 GeV have been produced in just a few centimeters [3H6] 
with 10 GeV stages being planned as modules of a high 
energy collider [7 , and detailed simulations are required 
to realize the promise of much shorter particle accelera- 
tors using this technique [8]. Such simulations challenge 
or exceed current capabilities, in particular for high en- 
ergy stages at GeV energies and beyond. 

The linear theory predicts that for the intense lasers 
(a>l) typically used for acceleration, the laser depletes 
its energy over approximately the same length L^ = 
A^/2Aq over which the particles dephase from the wake, 
where X p = ^nc 2 m/ e 2 n e is the plasma wavelength, Ao is 
the laser wavelength, c is the speed of light, and m, e and 
n e are respectively the electron mass, charge and density 
in the plasma pQ. As a result of beam dephasing and 
laser depletion, the maximum bunch energy gain scales 
approximately as the square of the plasma wavelength 
and the inverse of the plasma density, which implies that 
higher energy stages operate with longer plasmas, rend- 
ing computer simulations more challenging, as the ratio 
of longest to shortest spatial lengths of interest (plasma 
length/laser wavelength) rises. As a matter of fact, direct 
explicit multi-dimensional simulations of 10 GeV stages, 
which will operate in m-scale plasmas at order 10 1T /cc 
densities, have been considered until recently beyond the 
current state of the art [8j [9] . 

Recently, first principles Part icle-In- Cell modeling 
of laser-plasma wakefield accelerators using a Lorentz 
boosted frame of reference [10] have been shown to be- 
ing sped up by up-to three orders of magnitude in the 
calculations of stages in the 100 MeV-10 GeV energy 
range [8j [TTHT8] . Maximum obtainable speedups calcu- 
lated using linear theory predict that higher speedups are 
attainable, in the range of 4-6 orders of magnitude for 



stages in the energy range of 10 GeV-1 TeV respectively 
[19] [20]. Practical limitations have prevented reaching 
these speedups, including a violent high frequency nu- 
merical instability, limiting the Lorentz boost 7 below 

100 m nai ng mj [20] . 

We report for the first time direct explicit simulations 
of stages in the range of 0.1 GeV-1 TeV, using a Lorentz 
boosted calculation frame with gamma as high as 1300, 
verifying the performance and energy gain scaling [9j [21] 
of plasma accelerator stage with deep laser depletion into 
the 1 TeV range and providing the tools for detailed de- 
signs for upcoming 10 GeV experiments such as BELLA 
[22] . As we have shown in [10], the speedup provided 
by computing using a Lorentz boosted frame comes from 
the properties of space and time contraction and dila- 
tion of the Lorentz transformation. In this paper, the 
property of rotation of space-time of the Lorentz trans- 
formation is utilized to overcome the numerical instabil- 
ity that has arisen for boost values needed for reaching 
the maximal theoretical speedup. In conjunction with 
the development of novel numerical techniques that are 
described elsewhere p~9j [20] , this allows the simulations 
to approach the theoretically calculated speedups of 4-6 
orders of magnitude for 10 GeV-1 TeV stages, which in 
turn allows simulations of high energy plasma accelera- 
tors. 

Effect of the hyperbolic rotation in Minkowski space. — 
The effects of the Lorentz transformation on the laser 
and wake propagation through a 100 MeV laser plasma 
acceleration stage [9j[2T] is illustrated in space in Figure[l] 
and in space-time in Figure [2] taken from simulations 
using the Particle-In-Cell code Warp [23]. The Lorentz 
transformation can be described as a hyperbolic rotation 
in Minkowski space and its rotational effect is explicitly 
visible in Figure [2] 

Figure [l] shows surface renderings of the transverse and 
longitudinal electric fields respectively, as the beam en- 
ters its early stage of acceleration by the plasma wake, 
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FIG. 1: (color online) Colored surface rendering of the transverse (laser) and longitudinal (wake) electric fields from a 2-1/2D 
Warp simulation of a laser wakefield acceleration stage in the laboratory frame (left) and a boosted frame at 7 = 13 « jwake 
(right), with the beam (white) in its early phase of acceleration. The laser and the beam are propagating from left to right. 



from calculations in the laboratory frame and a Lorentz 
boosted frame at 7 = 13 (approximately the laser group 
velocity in the plasma column j g w 13.2). The two snap- 
shots offer strikingly different views of the same physi- 
cal processes: in the laboratory frame, the wake is fully 
formed before the beam undergoes any significant accel- 
eration, the laser (Fig. [I] top-left) is easily recognizable 
(i.e. its shape is only slightly distorted by the plasma) 
and leaves a visible imprint on the wake (longitudinal) 
field (Fig. [I]bottom-left); in the boosted frame, the beam 
is accelerated as the plasma wake develops, the laser (Fig. 
[l] top-right) is not easily recognizable (i.e. its shape is 
highly distorted by the plasma) and no evident imprint 
is left on the wake field (Fig. [I] bottom-right). 

The physics underlying the differences between 7 = 1 
and 7 = 13 views of the wake is illustrated by histories 
of the (transverse) laser field on the longitudinal axis, 
reported in Figure [2] from simulations performed using 
the laboratory frame and boosted frames at 7 = 5 and 
13. Simulations with 7 = 1 and 5 used a moving win- 
dow propagating at the group velocity v g of the laser 
in the plasma and the data are plotted in the frame 
of the (galilean) moving window. The simulation with 
7 = 13 did not use a moving window as it was per- 
formed very near the group velocity of the laser in the 
plasma j g ~ 13.2. The data from the boosted frame 
simulation at 7 = 13 is presented in the boosted frame 
as well as in the laboratory frame moving window (after 
Lorentz transformation) , allowing direct comparison with 
the simulation calculating with the laboratory frame. 
The calculation with the boosted frame at 7 = 13 was 
approximately 200 times faster than the calculation with 
the laboratory frame, as expected [10 . The agreement 
between the two is nonetheless excellent (comparing top- 
left plot to bottom-right plot in Figure [2|, confirming the 
accuracy of the calculation in the boosted frame. The 
group velocity of the wake in the plasma is always below 
the speed of light in vacuum while the phase velocity is 
always above it, resulting in oblique stripes in the labo- 



ratory frame plot. As 7 boost rises, the stripes rotate ac- 
cording to the rules of the Lorentz transformation, even- 
tually becoming nearly perpendicular to the time axis 
as 7 boost nears 7^ (bottom- left plot in Figure [2]), the 
laser group velocity approaches zero and the phase ve- 
locity approaches infinity. In effect, the laser oscillations 
that appear in the laboratory as spatial oscillations prop- 
agating in the plasma are transformed into time beating 
of the field for calculations in frames whose boost nears 
the laser group velocity. As discussed below, this effect 
has important consequences for the modeling of full scale 
stages at 10 GeV or above. 
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FIG. 2: (color online) Color plots of the laser field history 
on axis in a moving window propagating at the laser group 
velocity from simulations: (top-left) in the laboratory frame; 
(top-right) in a boosted frame at 7 = 5; (bottom- left) in a 
boosted frame at 7 = 13; (bottom-right) in a boosted frame 
at 7 = 13 plotted in laboratory coordinates after Lorentz 
transformation of the data. 
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Mitigation of a numerical instability. — Several numer- 
ical limits have restricted the boost performance in past 
simulations: laser initialization, statistics and a short 
wavelength instability. The first two are discussed else- 
where [19 and we concentrate here on the latter. 

A violent high-frequency numerical instability devel- 
oping at the front of the plasma column for boosts at 
7 > 100 in 2D and 7 > 50 in 3D was reported by various 
authors p~5j [16j [24]. The presence and growth rate of 
the instability was observed to be very sensitive to the 
resolution (slower growth rate at higher resolution), to 
the amount of damping of high frequencies, to smooth- 
ing of short wavelengths, and to the boost value (stronger 
instability at higher boost). An extensive set of testing 
was performed with Warp to investigate simulations of 
downscaled 100 MeV and full scale 10 GeV LPA stages 
[20] . One of the key findings was that frames with higher 
boosts (up to 7 P ) allow for higher levels of filtering and 
damping than is possible in other frames for the same 
accuracy, allowing mitigation of the instability. This is a 
direct consequence, and benefit, of the hyperbolic rota- 
tion effect from the Lorentz boost seen in Figure [2] 
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FIG. 3: (color online) Spectral content history of the laser 
field on axis as it propagates through the plasma column for 
100 MeV (top) and 10 GeV (bottom) stages, in the laboratory 
frame (left) and near the laser group velocity frames at 7 = 13 
(top-right) and 130 (bottom-right). 



in a narrow band around the nominal laser wavelength, 
then spreads due to dispersion and depletion effects as 
the laser propagates through the plasma. In the frame 
of the laser group velocity, much of the spectral content 
is localized initially at wavelengths that are several times 
the nominal laser wavelength (in vacuum), then progres- 
sively fills lower portions of the spectrum. As in practice 
the longitudinal numerical resolution is set relative to 
the vacuum laser nominal wavelength, a higher level of 
filtering (or damping) is acceptable for simulations in a 
boosted frame than for those in the laboratory frame. 
Furthermore, the higher the boost, the longer the wave- 
lengths with substantial spectral content in the early part 
of the laser propagation (comparing spectral content be- 
low diagonal in plots from right column of Figure [3|, 
meaning that higher levels of filtering are allowable at 
higher boosts were the growth of the numerical instabil- 
ity is the strongest. 
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FIG. 4: (color online) Verification of the scaling [9j [21] of 
electron beam energy gain versus longitudinal position (in 
the laboratory frame) from direct simulations at n e — 10 19 cc 
down to 10 15 cc (energy gains from 0.1 GeV to 1 TeV), using 
Lorentz boosted frames of reference at 7 between 12 and 1300, 
in 2-1/2D (top) and 3D (bottom). Energies from the lowest 
energy stages were scaled to the highest energy stage to allow 
for direct comparison. 



The spectral content history of the laser field on axis as 
it propagates through the plasma column is given in Fig- 
ure [3] for 100 MeV and 10 GeV stages, in the laboratory 
frame and in the laser group velocity frames at 7 = 13 
and 130 respectively. Laser depletion occurs at times 
Td = Ld/c at respectively & 3.1 ns and ~ 3.1 ps 
for 100 MeV and 10 GeV stages. In the laboratory frame, 
the spectral content of the laser is concentrated initially 



Modeling of up to 1 TeV stages. — Mitigation of the 
instability allowed simulations of stages with gain energy 
as high as 1 TeV in 2-1/2D and 100 GeV in 3-D, us- 
ing Lorentz boosted frame with 7 as high as 1,300 (see 
Figure [4]), offering for the first time direct verification of 
the scaling of plasma accelerators into the 1 TeV range 
for deeply depleted stages. Simulations with the high- 
est boost necessitated filtering of short wavelength over 
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a wider band, and in agreement with the observations of 
the previous section, the accuracy was not compromised. 
The highest level of smoothing was needed for the 1 TeV 
case, explaining the deviation past 1 km. This deviation 
is of little importance in practice, where one is mostly 
interested in the beam evolution up-to the peak energy 
point. The differences at n e = 10 19 cc can be attributed 
to the effects from having only a few laser oscillations 
per pulse. The theoretical speedup [19] of the full scale 
100 GeV class run is estimated to be over 100,000. As- 
suming the use of a few thousands of CPUs, a simulation 
that would have required an impractical several decades 
to complete using the laboratory frame, was completed 
in only four hours using 2016 CPUs of the Cray system 
at NERSC. The speedup of the 2-1/2D 1 TeV stage is 
estimated to be over a million. 

The boosted frame Particle-In-Cell technique accu- 
rately resolves the wavelength shifting and broadening 
that occurs as the laser depletes, offering advantages over 
other models (for example envelope, quasistatic) while 
providing the speed required for direct simulation of 10 
GeV and beyond laser plasma accelerators to accurately 
model laser and beam transverse oscillations. It is being 
applied to the direct simulation of 10 GeV beam loaded 
stages for detailed designs of experiments on new lasers 
such as BELLA [22] (see Figure |5|, as well as next gener- 
ation controlled laser plasma accelerator stages and col- 
lider modules [7]. 




FIG. 5: (color online) Snapshot from a 10 GeV LPA stage 
boosted frame simulation. The image shows an externally 
injected electron bunch (rainbow color) riding a density wake 
(blue) excited by an intense laser pulse (red), propagating in 
a 0.65 m long plasma channel. The laser pulse ( 40 J in 67 fs), 
focused to 90 /xm spot size at the entrance of the channel, has 
reached the end of the plasma channel. The electron bunch 
energy has reached up to ~10 GeV. 

In summary, direct simulations of stages in the range 
of 0.1 GeV- 1 TeV have been performed using the Lorentz 
boosted frame technique, verifying for the first time the 
performance of plasma accelerators into the 1 TeV range 
for deeply depleted stages. This has been possible thanks 
to effects of the hyperbolic rotation in Minkowski space 
of the laser propagation in the plasma column which has 
been key, in conjunction with the development of novel 
numerical techniques (described elsewhere), in allowing 



successful mitigation of a violent numerical instability 
that was limiting the boost performance in past simu- 
lations. As a result, the maximum theoretical speedup 
of over a million for a 1 TeV stage was realized, which 
is three orders of magnitude higher than was possible 
previously. The new developments offer unique highly 
efncent tools for the detailed designs of experiments on 
new lasers such as BELLA [22] . 
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